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Abstract. Partial differential equations with discrete (concentrated) state-dependent 
delays are studied. The existence and uniqueness of solutions with initial data from a wider 
linear space is proven first and then a subset of the space of continuously differentiable 
(with respect to an appropriate norm) functions is used to construct a dynamical system. 
This subset is an analogue of the solution manifold proposed for ordinary equations in 
^ ■ [H.-O. Walther, The solution manifold and C 1 -smoothness for differential equations with 
state-dependent delay, J. Differential Equations, 195(1), (2003) 46-65]. The existence of 
a compact global attractor is proven. 



1 Introduction 

The partial differential equations (PDEs) with delays have attracted a lot of attention dur- 
ing the last decades as many processes of the real world (like an automatically controlled 
furnace, bi-directional associative memory (BAM) neural networks, reaction-diffusion pro- 
cesses) can be described by such kind of equations. Studying these equations is based on 
the well-developed approaches to the ordinary differential equations (ODEs) with delays 
[THE,!!] and PDEs without delays [Stl5l H5lfllj . Under certain assumptions both types of 
equations describe a kind of dynamical systems that are infinite-dimensional, see [2j [301 [6] 
and references therein; see also [311 [U El [3] and to the monograph [37] that are very close 
to this work. 
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In many evolution systems arising in applications the presented delays are frequently 
state- dependent (SDDs). The theory of such equations, especially the ODEs, is rapidly 
developping and many deep results have been obtained up to now (see e.g. [32j [33j [34j 
fT6l [18] 135] and also the survey paper [12] for details and references). The underlying 
main mathematical difficulty of the theory of PDEs with SDDs lies in the fact that the 
functions describing state-dependent delays are not Lipschitz continuous on the space of 
continuous functions - the main space, on which the classical theory of equations with 
delays is developed. This implies that the corresponding initial value problem (IVP) is 
not generally well-posed in the sense of J. Hadamard [HJ [9]. 

The partial differential equations with state-dependent delays were first studied in 
[2T] (the case of distributed delays, weak solutions), [13] (mild solutions, infinite discrete 
delay), and |22j (weak solutions, finite discrete and distributed delays). An alternative 
approach to the PDEs with discrete SDDs is proposed in |24j . 

This paper is a continuation of the work |25j and its goal is to study the approach 
used for ODEs with SDDs [32], [33], [12] in the case of PDEs. The main idea lies in finding 
a wider space Y D X such that a solution u : [a, b] — > Y be a Lipschitz function (with 
respect to a weaker norm of Y), and constructing a dynamical system on a subset of the 
space C([a,b];Y). It should be emphasized that the dynamical system is constructed on 
a metric space that is nonlinear. More precisely, the existence and uniqueness of solutions 
with initial data from a wider linear space is proven first and then a subset of the space 
of continuously differentiable (with respect to an appropriate norm) functions is used 
to construct the aforementioned dynamical system. This subset is an analogue of the 
solution manifold proposed in [33], see also [12] . We use the same class of non-local in 
space variables nonlinear PDEs as in [25J. 

The paper is organized as follows. The section 2 is devoted to the formulation of the 
model. The proof of the existence and uniqueness of (strong) solutions for initial functions 
from a Banach space forms a main part of the section 3. In the section 4, an evolution 
operator S t is constructed and its asymptotic properties in different functional spaces are 
investigated. The dissipativeness is obtained in a Banach space, while the existence of a 
global attractor is proven on a smaller metric space (the solution manifold). The choice 
of this smaller space is different from that proposed in [25] . 
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2 The model with discrete state-dependent delay and 
preliminaries 



Consider the following non-local partial differential equation with a discrete state- 
dependent delay 77 

d 

—u(t,x) + Au(t,x) + du(t,x) =b([Bu(t - T](u t ), -)}(x)) = (Fi(w t ))(x), x eQ, (1) 

where A is a densely-defined self-adjoint positive linear operator with domain D(A) C 
L 2 (Q) and compact resolvent, which means that A : D(A) — » L 2 (f2) generates an analytic 
semigroup, f2 C R no is a smooth bounded domain, B : L 2 (fi) — > L 2 (Q) denotes a bounded 
operator that will be defined later, b : R — ?• R stands for a locally Lipschitz map, d G 
R, <i > 0, and the function 77 : C([— r, 0]; L 2 (fi)) — >■ [0,r] C R+ denotes a state- dependent 
discrete delay. Let C = C([— r, 0]; L 2 (fi)). Norms defined on L 2 (Q) and C are denoted by 
|| • || and || • ||c, respectively, and (•, •) stands for the inner product in L 2 (Q). As usually, 
u t = u t (9) = u(t + 6) for 9 G [-r, 0]. 

Remark 1. The operator i? may for example be of the following forms (linear operators) 



where / : Q — > R is a smooth function and I G C^°(fl). In the last case the nonlinear 
term in ([!]) is of the form 




(2) 



or even simpler 




(3) 



(Fi(ut)) (x) = b ( [ u(t - v(ut),y)f(x - y)£(y)dy 



(4) 



□ 



Consider the equation (CQ) with the initial condition 



U|[-r,0] = <P 



(5) 



and let 



H ee |v9 G C7([-r,0];£)(A-3)) | ^(0) G 



(6) 
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Let further 



\(P\\h= max ||A~3yj(s)|| + ||Ai^(0)|| 

se[-r,0] 



be a norm defined on the space H and D(A a ) denote the domain of the operator A a . In 
the sequal the following assumptions will play an important role. 

(HI. 77) The discrete delay function r] : H — > [0, r] is such that 

3L V > 0, 3q > such that VV>, ip G H =>• 

\n{<p) - vw\ < l v (quarto) - + f u-Hm - mw^j 2 (7) 

(H.B) The following Lipschitz property of the operator B holds. 

3L B > such that Vu,v G D(A~?) => \\Bu - Bv\\ < L B \\A~^(u - v)\\ (8) 



Remark 2. Under the assumption that for all (almost all) x G Q =>- /(■ — x)£(-) G D(A2 
and u G D(A~^), the term of the form ([3]) implies that 

| <«,/(• -ar)*(-)>| < \\A-5u\\\\A3f(.-x)£(-)\\, 

which gives 

/ u(y)f(y-x)£(y)dy\ 2 dxY < \\A^u\\ ( [ \\A^f(- - x)£(-)\\ 2 dx 



Hence, the property (H.B) (see ©) holds with L B = (/ Q \\A*f(- - x)£(-)\\ 2 dxj * The 

same arguments hold (with L B = ^ J n \\A^f(x, -)\\ 2 dx^j 2 ) for a more general term of the 
form (ED. □ 



Let now the following space 

i\\A-U^)-r(t))\[ 



C={ V Z C([-r, 0];D(A-i)) | sup { <P<t))\\ ] < +00; m e D{M) j> , 



' (9) 

with the natural norm 

II II - \\A~h ( Ml , | P~^( g )-l(X>)ll l . / n xii , ln x 

MLc = max p 2^(s)||+sup<^ , ) + \\A*<p(p)\\ 10 

se[-r,0] s ^ I |S - t| J 



be defined. For any segment [a,b] C R (c.f. (J5])) and any Lipschitz-on-[a, b] function ip, 
let 

||M|| M1 = sup j H- 4 "*'^ ~ P(t))l1 : t; s, t e [a,b)\ (11) 
denote its Lipschitz constant and let |||<£>||| = 1 1 \(p\ \ \ [- r ,o]- Then the following lemma holds. 



Lemma 1. Let the assumptions (HI. rj) and (H.B) hold (see ©, (|HJ)) and let the function 
b : R — > R is Lipschitz and bounded (|6(s)| < Mb for all s E R). Then any two functions 
(p E £,ip E H (with H and £ defined in (jB]) and Q) the nonlinear ity F satisfies 



(gp 1/2 (^(0) - #)))|| + \\A- l '\ v - i>)\\ c ) , (12) 



where 

L Fl [£] = LfeL^v^max {l; max{l; a/t}} (13) 
and L Fl [£] is used in (fT2|) with 

, = s IIMH s sup { "■ 4 " /2 <^)-^) ) " : t; s,t e [-r,0]} . 

Proof of Lemma 1. Using the Lipschitz property of b and S (see (H.B)), it follows 
that 

\\F x (<p) - F x ty)\\* = f IKIB^-^)^)) - b([B^}(- V ^),x))\ 2 dx < 
Jn 

<L 2 f |[^^](-r 7 (<^),a;)-[^](-r 7 (V), a:)| 2 ^ = 1 1 [^^] ( -77 (^) , -) - [^^] (-^(V) , OlP < 
Jn 

< L 2 b L%\\A-y* M-rtvO) - ^(-r/W) ± vi-vim II 2 < 

< 2L 2 L| (HA- 1 / 2 foOKvO) - ^W)} || 2 + \\A- 1 ' 2 (<p - VO H 2 ,) • 
Next, tp E C implies that there exists L v = |||</?||| > 0, (see (TTTJT) . (TTTT) ) such that 

WA-^ipis 1 ) - V (s 2 ))\\ < L^s 1 - s 2 \, VsVe[-r,0]. (14) 
Hence, (fT4"]l and (HI. 77) give 

iii^-i^voir^ 
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< 



< 2l 2 l| 



2 c d6)^ 



<2LlLl\LlLl( q \\A l l\m-m)\\ 2 + J \\A-i{<p{0)-il>{0))\\ a M) + 

+ \ \A~* (ip- | 
%A (q \\A^(0) - + r \\A-H<p 

+ \ \A~% 

< 2L\L\ max {1; L%L\ max{l; r} } [q \\A^ 2 {ip(0) - ^(0)) 1 1 2 + \\A~i {<p - ^ 

The last estimate and using the formulas ^/max{|a|; \b\} = max{y|a|; y/\b\} and 
\/a 2 + b 2 } < \a\ + \b\ give (fT2l . (Fl3|) . which completes the proof. □ 



< 



3 The existence and uniqueness of solutions 



As in [25] we need the following 



Definition 1. A vector-function u(t) G r, T]; D(A^ 1 ^ 2 )) n C([0, T}; D(A 1 ^ 2 )) n 

L 2 (0,T;L>(A)) with derivative u(t) G L°°(0,T; D(A- 1/2 )) is a (strong) solution to the 
problem defined by (OQ) and (jSJ) on [0, T] if 

(a) u(0) = <p(6) for G [-r,0]; 

(b) Vv G L 2 (0, T; L 2 (Vl)) such that u G L 2 (0, T; D(A^ 1 )) and t»(T) = => 



rp rp 

(u(t),v(t))dt + [ {A 1 ' 2 u{t),A 1 l 2 v{t))dt = 
i) Jo 

= (<p(P),v(0))+ f (F^-d-u^^dt. (15) 
Jo 



Now we prove the following theorem on the existence and uniqueness of solutions. 



Theorem 1. Let the assumptions (HI.?]) and (H.B) hold and let the function b : R — > R 
be Lipschitz and bounded, i.e. \b(s)\ < M b for all s G R. Let further <p G C be a given 
initial condition. Then the problem defined by ([1]) and ([5]) has a unique solution on any 
time interval [0,T] such that tt G L 2 (0, T; L 2 (fi)). 
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Remark 3. Notice that ip does not assume ip G L 2 ([—r, 0]; D(A)). However, the definition 
of a strong solution implies that 



Ut eL 2 ([-r,0];D(A)), Vt > r. 



(16) 



□ 



Proof of Theorem 1. We follow the proof of Theorem 1 in [25]. Notice that the as- 
sumption (HI. 77) is slightly more general than the assumption (H.77) in [25]. This implies 
some changes in the proof of the uniqueness of solutions. 

Let {e/ c }^ 1 denote an orthonormal basis of L 2 (Q) such that Ae k = X^k, < Ai < 
... < Afc — > +00 and consider the Galerkin approximate solution u m = u m (t,x) = 
Y^k=i9k,m{t)^h of order m such that 



Wk = 1, . . . ,m, g k)7n G C 1 (0, T; R) fl L 2 (— r, T; R) with gk,mif) absolutely continuous. 

The system (II7|) is a system of (ordinary) differential equations in R m with a con- 
centrated (discrete) state-dependent delay for the unknown vector function U(t) = 
(gi,m(t), . . . ,g m , m {t)) (for the corresponding theory see [331 El] and also a recent review 



The key difference between equations with state-dependent and state-independent 
(concentrated) delays is that the first type of equations is not well-posed in the space 
of continuous (initial) functions. To get a well-posed initial value problem, it is better 
[33| 13"^ IT2] to use a smaller space of Lipschitz continuous functions or even a smaller 
subspace of C l ([-r, 0];R m ). 

The condition (p G C implies that the function t/(-)|[- r ,o] = Pm<-p{-), which defines 
initial data, is Lipschitz continuous as a function from [— r, 0] to R m . Here P m is the 
orthogonal projection onto the subspace span {ei, . . . , e m } C L 2 (Q). Hence, we can apply 
the theory of ODEs with discrete state-dependent delay (see e.g. [12]) to get the local 
existence and uniqueness of solutions to ffT7|) . 

Next, we will get an a priory estimate to prove the continuation of solutions u m to ( 117]) 
on any time interval [0,T] and then use it for the proof (by the method of compactness, 
see [15]) of the existence of strong solutions to <Q and (jS]). To that end, multiply the first 



(ii m + Au m + du m - Fi«), e k ) = 0, 
(u m (6),e k ) = (<p(9),e k ), V0e[-r,O] 



(17) 



H21). 
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equation in (fT7j) by \kgk,m and sum for k = 1, . . . , m to get 

l^pl/2 u m W ||2 + p u m (t) ||2 + d . ||^l/2 um(t) ||2 = (p^f), A M ™(t)> < 

<i||P m F«)|| 2 + ^||^ m (t)|| 2 . 
As the function b is bounded, H-F^w™)!! 2 < M 2 |f2| (here \Q\ = J Q Ida;), which gives 

^\\A 1/2 u m (t)\\ 2 + ||A M m (t)|| 2 < M 6 2 |0|. (18) 

Integrating (ITS]) with respect to i and using the relationships <p(0) G D(A 1/2 ), w m (0) = 
P m y>(0) G D(v4 1/2 ), \\A l ' 2 u m {<d)\\ = 1 1^/^^(0)1 1 < IIA^VCO)!!, we get an a priory 
estimate 

\\A 1/2 u m {t)\\ 2 + f \\Au m (r)\\ 2 dr <\\A 1/2 if{0)\\ 2 + M 2 \Q\T, Vm,VtG[0,T]. (19) 
Jo 

The above relationship (|T9|) means that 

{u m }™ =1 is a bounded set in L°°(0, T; D(A 1 ' 2 )) n L 2 (0, T; D(A)). 
Using this fact and ( fTTj) . it follows that 

{u m }™ =1 is a bounded set in L°°(0, T; £>(A- 1/2 )) n L 2 (0, T; L 2 (fi)). 
Hence, the family {(ii m ; u m )}™ =1 is a bounded set in 

Z x = (L°°(0, T; D(A 1/2 )) n L 2 (0, T; x 

x (L°°(0, T; ^(A- 1 / 2 )) n L 2 (0, T; L 2 (fi))) . (20) 

Therefore, there exist a subsequence {{u k ;u k )} and an element (u;u) G Zi such that 

{(u k ;u k )} *-weak converges to (u;u) in Z\. (21) 

The proof that any *-weak limit is a strong solution is standard. To prove the property 
u(t) G C([0,T];D(A 1 / 2 )), we use the well-known (see also [14, thm. 1.3.1]) 

Proposition 1. (Proposition 1.2 in [26]). Let V denote a dense Banach space that is 
continuously embedded in a Hilbert space X and let X = X* so that V > X c -» V*. Then 
the Banach space W p (0, T) = {u G L p (0, T;F):«6 L 9 (0, T; 1/*)} (here p" 1 + g" 1 = 1) is 
contained in C([0, T];X). 
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In our case X = D(A l l 2 ), V = D(A), V* = L 2 (Q), p = q = 1/2. 

Now we prove the uniqueness of solutions. Using the fact that if G C, the definition 1 
of a solution v, and v(t) G L°°(0, T; D(A^ 1 ^ 2 )) (see (|2TT) ) . it follows that for any such a 
solution v and any T > 0, there exists L Vi t > 0, such that 

IIA-^Ms 1 ) -v(s 2 ))\\ < L VtT \s l - s\ Ws\s 2 G [-r,T]. (22) 

In the light of ( ITTT) . let L V) t = \ \\v\\\i- r ,T]- 

Consider any two solutions u and u of fl}, (JSJ) (not necessarily with the same 
initial function). The standard variation-of-constants formula u(t) = e~ At u(0) + 
J*e- A ^F(u T )dT and the estimate p a e" M || < (f) a e~ a (see e.g. [gj (1.17), p.84]) 
give 

\\A^ 2 (u(t)-v(t))\\ <e- Xlt \\A^ 2 (u(0)-v(0))\\+ f \ \A^ 2 e~ A ^ \ \ \\F{u t )-F{v t )\\cIt < 

Jo 

< e-^\\A l / 2 (u(0)-v(0))\\ + J* (^y 2 e- 1/2 \\F(u T )-F(v T )\\dr } (23) 

1 /2 

as \\A l / 2 e~ A ^\\ < (0} e~ 1 / 2 and similarly, 

\\A- l ' 2 {u t -v t )\\ c < \\A- l ' 2 (u, -v )\\c+ f\\F{u T ) - F(v T )\\dr. 

Jo 

The last estimate and (1231) give (just the case when q = 1 is shown for the purpose of 
clarity) 

\\A^ 2 {u{t) - v{t))\\ + \\A- l ' 2 (u t - v t )\\c < e- Xlt \\A^ 2 (u(0) - v(0))\\ + 

+ \\A- 1 / 2 (u -v )\\c+ f {l + {2e{t-T))- l l 2 }\\F{u T )-F{v T )\\dr. (24) 

Jo 

It follows, from Lemma [IJ that 

\\F{u t )-F{v t )\\ < L Fli%T (qWAWW) -v(t))\\ + \ \A- X ' 2 (u t - v t )\\ c ) , (25) 

where L Fl<VtT is defined in the same way as L Fl in (fl3"|) . just with £ = L v T instead of L v 
- see (USD and (T22l). 

= LfoLsV^max {l; L Vj tL v max{l; a/t}} . (26) 



J Fi,v,T 



= L 



Ft 



J v,T 



9 



It should be emphasized how the Lipschitz constant L v ^ = \\\v\\\\^ r T\ of a strong 
solution v is taken into account in (126]) (see (1221) and dill)). 



Let 

0(t) = ||i4 1/2 («(t) - v(t))\\ + WA- 1 / 2 ^ - v t )\\ c , 
Then the relationships (124p and ( 1251) lead to the following estimate 

g(t) < g(0) + / {1 + (2e(t - r))- 1 / 2 } ■ #(t) dr 

Jo 

Lemma 2 (Gronwall). Let u, a G C[a, > 0,/3 is integrable on [a, 6] and 

u(t) < a(t) + I 0(t)u(t) dr, a<t<b 

J a 

Then ^ ^ 

u{t)<a{t) + J /3(r)a(r)exp|^ dr, a<t<b 

Moreover, if a is non- decreasing, then 

u(t) < a(t)exp|y /3(s)ds|, a < t < b. 

It follows, from the above lemma and equality J*„ (t — r) _1 / 2 <ir = 2t 1 / 2 , that 

g(t)<g(0) exp^L FltVtT J {l + (2e(t - s))" 1/2 } dsj < 

< gf(0) exp |l Fi)WiT + 
which implies, Wt e [0,T], that 

H^Kt) - «(f))|| + HA- 1 ^^ _ Vt) \\ c < 

< (HA^KO) - «(0))|| + HA- 1 ^^ _ Vo) \\ c ) , (27) 



where 



^i,d,t = exp <( Lir lit))T • ( T + J ^ ) } , (2£ 



see (1261) for the definition of Lf 1>v> t = L Fl [L V) t\- This proves the uniqueness of the 
solution to (pQ) and (JSJ) , and completes the proof of theorem 1 . □ 
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4 Asymptotic properties of solutions 



This section is devoted to studies of the asymptotic behavior of solutions in different 
functional spaces. We define first (in a standard way) the evolution semigroup St : C — > C 
(the space £ is defined in fl9])) by the formula 

Sup = u t , t> 0, (29) 

where u(t) is a unique solution to the problem p]) and fl5]) (see definition 1). 

The estimate (f27|) means the continuity of the evolution operator St in the norm of 
the space H (see fl6])), i.e. 

\\S t <p-S t ip\\H<E Fl , v , T \\<p-ip\\H for all t€ [0,T]. (30) 

The aim now is to get a more precise estimate, e.g. the continuity of St in the norm of 
the space £ (see (jHJ), (jTU]) ). Consider the definition of the Galerkin approximate solution 
(see ( JIT]) ). It gives 

||A- 1/2 (M m (t) - i> m (t))|| < ||A 1/2 ( M m (t) - v m (t))|| + d\\A~ 1/2 (u m (t)- 

-v m (t))\\ + \mo-FiW\\ 

and Lemma [1] implies 

1 1 A- x / 2 (u m (*) -v m (t))\\ <(l + d + L F1 ){\\A 1/2 (u m (t) -v m (t))\\+ 

+ \\A-V 2 (u?-v?)\\ c }. 

An analogous estimate for a solution to the problem (JTJ and (jSJ), can be obtained from 
(1211) and the following 

Proposition 2. [38| Theorem 9] Let X be a Banach space. Then any *-weak conver- 
gent sequence {wj.}™ =1 G X* *-weak converges to an element Woo G X* and ||Woo||x < 
liminf^oo ||w n |U- 

More precisely, 

ess sup t€l0T] \\ A- 1 / 2 {ii{t) - v{t))\\ < (l + d + Lm) sup {HA 1 ' 2 ^) - v(t))\\+ 

te[o,r] 

+ ||A- 1 / 2 (« t -^)|| C7 } 
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The last estimate and relationship ( 12 7p imply 



^Pte[o,T]U- 1/2 (u(t)-v(t))\\ < 

< (l + d + L F1 )E Fl>v>T (\\A^ 2 (u(0)-v(0))\\ + \\A- l l\u Q -v )\\ c ) (31) 



CSS 



Hence, see (fTTj) . 

WW- v\\\ [m <{l + d + L Fl )E Fl ^ T (\\A^ 2 (u{0) - v(0))|| + WA-^iuo - v )\\ c ) 

From that and (1271) . it follows that 

\\ut-v t \\c<(2 + d + L F1 )E Fl:V:T \\uo-v \\c, VtG[0,T], (32) 

which finally means that for any T > there exists a constant C F > such that Vt G [0, T] 
it gives 

IK - v t\\c = \\S t <p - S t ip\\c < C T \\<p - ip\\c, \/(p,ipe£ (33) 

The last inequality means the continuity of the evolution operator St in the norm of the 
space £ (see (129]) and compare with (|30|) ). 

Remark 4. It should be noted that the evolution operator and, more generally, the time- 
shift is not a (strongly) continuous mapping in the norm of the space £ (see (jHJ)). This 
can be illustrated by the following simple (scalar) example. 
Consider the space 

Cip([-r,T\;R) = \ v : [-r,T] -»■ R : sup ( ^ y , , * ^ t;s,te [-r,T}\ < oo 



\8-t i 

and analogously define the space dp ([— r, 0]; R) with the natural norm 

IHU<p= max \v{6)\ + sup j^V — rr^> s ^ *5 s ^ e h r ' °1 1 • 

6e[-r,0] { \S — t\ J 

The (strong) continuity of the time-shift means that 

Vv G Cip([-r,T\;R) and Vt G [0,T] =J> lim | |u m - v t | | £ip = 0. (34) 

h— >0 

Obviously, when t = one considers h + , while for t = T, the case /i — > 0~ should be 
investigated. 

To prove the claim, we must show that (1341) does not hold, i.e. 
3v G £zp([— r, T; R) and 3t £ [0>^1 f° r which lim |K +/, — ^tollr™ 7^ 0- (35) 
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Thus, consider the case tn = 0, h — > + and the function 



v{t) = 

It can be seen that v to — v = and 

vt +h = v to+h (0) 



0, te[-r,0) 

t, te{o,T] 



0, G [-r, -/i] 

+ 5e(-/i,o] 



Hence, | \v to+h - v to \\ Cip = \ \v to+h \\ C i P = h + 1 and finally lim | \v to+h - v to \\a P = lim (h + 

h— >0+ h—>-0+ 

1) = 1 / 0, which means that (JM|) does not hold. In the space £, we would proceed 
analogously. □ 



Remark 5. In the same way as in the previous remark one can show that the time-shift is 
not a (strongly) continuous mapping in the topology of L°°(—r, 0). One could consider the 
|o, te[-r,0] 

function v (t) = < and to = to show that lim \\vh — Uo|U°°(-rO) = 1^0. 

[i, te(o,T] 

By the way, v = 4u, where, as usually, the time-derivative is understood in the sense of 
distributions. □ 

The above remarks show that despite of the existence and uniqueness of solutions in 
the space C and even strong continuity of the evolution operator St in the norm of £ (see 
( 133|) ). the pair (St; £) does not form a dynamical system since St is not strongly continuous 
as a mapping of time variable. 

The methods, developed for ordinary delay equations in [33] suggest to restrict our 
considerations to a smaller subset of the space of Lipschitz functions. In this paper we 
follow this suggestion and consider the evolution operator St on the following subset of C 

X = { (p G C\[-r, 0}; D(A~ 1/2 )) such that 

<p(0) G D(A 1/2 )) and 0(0) + A<p(0) + d<p(0) = F x (<p) } C C. (36) 
Here the equality 0(0) + A(p(0) +dip(0) = F±(ip) is understood as an equality in D(A~ 1 ^ 2 ). 



13 



Remark 6. The set X is an analogue of the solution manifold introduced in [33] for the 
case of ODEs with state-dependent delays. □ 



To show that the set X is invariant under the evolution operator St, we first have to 
establish an additional smoothness property of the solutions of problem (Op), p]). 

Lemma 3. For any (p E C\[-r, 0]; D(A~ 1/2 )) such that (p(0) E D(A l/2 )), the solution to 
CQ), © (which is given by TheoremQ]) has the property (c.f. Proposition [1] and Theorem[T]) 

u E C([0, T]; D(A- 1/2 )), VT > 0. (37) 

Remark 7. We do not assume cp(0) E D(A), just cp(0) E D^A 1 / 2 ), so we cannot directly 
use [TJH Theorem 3.5, p. 114]. 

Proof of Lemma 0. By Proposition 1 and Theorem 1, for any E C 1 ([—r, 0]; D(A~ l l 2 )) 
such that y?(0) E D(A 1 ^ 2 )), there exists a unique solution u(t) E C([—r, T]; D(A^ 1 ^ 2 )) n 
C([0, T]; ^(A 1 / 2 )). This property and Lemma [1] then imply the continuity of the function 

p(t) = F^ut) E C([0,T];L 2 (n)). (38) 

Consider the following auxiliary linear system without delay 

v(t)+Av(t) + dv(t)=p(t), t>0, 
v(0) = (p(0) E D(A X / 2 ) 

In the same way as in fTl7|) . the Galerkin approximate solution v m = v m (t,x) = 
Y^k=i9k,m(t)e.k of order m to ( |39l) can be defined such that 

(v m + Av m + dv m - p(t), e k ) = 0, t > 0, 
(t; m (0),e fc ) = (^(0),e fe ), Vk = l,...,m. 

where gk, m £ C 1 (0,T; R) fl L 2 (— r, T; R) and gk, m {t) is absolutely continuous. 

The difference between approximate solutions u m and i> m lies in that v m are solutions 
just to linear system f T40|) . So, for any two approximate solutions v n and f m (solutions to 
( 1401) of different orders n and m), one has gk,n(t) = gk, m (t), which is denoted by (?fc(t). 

Multiply ( l4"0l) by A^^fc and sum for fc = n + l,...,n + p (pis any positive integer) to 

get 

(v n+p {t) - v n {t), A{v n+P {t) - v n {t))) + \ \A(v n+p (t)-v n (t))\\ 2 + 

+ d (v n+p (t) - v n (t), A(v n+P (t) - v n (t))) = ((P n+P - P n )p(t), A(v n+P (t) - v n (t))) 
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It should be recalled that, see the proof of Theorem [Q, P m is the orthogonal projection 
onto the subspace span {ex, . . . , e m } C L 2 (fl). Hence, 



2dr 



< \\(P n+p -Pn)p(t)\\ ■ \\A{v"**{t)-V»{t))\\ < l\\(Pn +p -Pn)p(t)\\ 2 + 

+ \\\A(v n+p (t)-v n (t))\\ 2 



which gives 

— II 
eft" 



A h(v n+ r(t) - v n (t))\\ 2 + - w B (t))H 2 < \\(P n+P ~ PnW)\\ 2 . 

Integrating the last estimate results (Wt G [0,T]) in 

\\A^{v n+p {t) -v n {t))\\ 2 + [ \\A{v n+p {r) - v n {T))\\ 2 dr < 

Jo 

< ||^(^(0) -«»(0)) 1 1 2 + f\\(P n+P 'Pn)p(r)\\ 2 dr< 

Jo 

<\\(P n+p -P n )A^(0))\\ 2 + [ \\(P n+p - P n )p(r)\\ 2 dr. 

Jo 

Summing up, the above estimate, the fact that (p(0) G D(A 1 ^ 2 )), the strong conver- 
gence \\I — P n \ \ — > for n — > oo, and ( |38l) imply that 



the sequence {v n }™ =l is a Cauchy sequence in C([0, T]; D(A? )). (41) 

Now our goal is to show that the sequence {v n }^ =1 is a Cauchy sequence in 

C([0,T];D(A-V 2 )). So, multiply first (SOD by to get X^g k (t) = -\\g k {t) - 

dX~Jgk(t) + (\f 2 P (t),e k ). This gives \^(g k (t)) 2 < 3\ k {g k {t)) 2 + 3d 2 {g k {t)f + 
_ i 

3\(\ k 2 p(t), e k )\ 2 . The sum for k = n + 1, n + p reads 

\\A-^{v n+p (t) - v n {t))\\ 2 < 3\\A^(v n+p (t) -v n (t))\\ 2 + 

+ 3d 2 1 1 {v n+p {t) - v n {t)) 1 1 2 + J- 1 1 {P n+P - P n )p(t) || 2 < 

A n+l 

^ 3 ( 1 + il") M V +P W - ^ n W)H 2 + ir-H J - ^H 2 ib(*)H 2 - 

V A n+lJ A n+l 

The last estimation together with ( 14~TT) give that 

the sequence {i; n }^l 1 is a Cauchy sequence in C([0, T]; D(A~ 5 )). (42) 

Thus, there exists a unique solution v(t) {v = lim^oo v n ) to the linear system ( |39l) . which 
satisfies v G C([0, T]; D(Aa)) and t> G C([0, T]; D(A~3)). 
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On the other hand, the nonlinear delay system (fTl) , (jSJ) with the initial function p> has 
also a unique solution. From the construction of p(t) (see fl38|) ). it follows that u(t) = v(t) 
for all t G [0,T], which gives (j37l) and completes the proof of Lemma [3l □ 

Lemma [3] particularly shows that the set X, defined by (|36|) . is invariant under the 
evolution operator S t (see ( 129]) ). This fact allows to define an evolution operator (denoted 
again by S t ) S t '■ X — > X in the same way as in (129|) . Now, if the natural norm 

M\x= max ||A- 1 /V(s)||+ max p" 1 / 2 ^)!! + II^VWII 

se[-r,0] se[-r,0] 

on X is taken into account, then Theorem [TJ Lemma [31 and Proposition [TJ give the 
continuity of S t with respect to t in the norm of X. Hence, (S t ; X) defines a dynamical 
system. 

Now we will pay attention to the long-time asymptotic behavior of the constructed 
evolution semigroup S t : X — > X. 

Theorem 2. Using the above notation and under the assumptions of Theorem [H 
the dynamical system (St,X) is dissipative. If, in addition, q = in (HI.?]), then 
(S t ,X) possesses a compact global attractor A, which is a bounded set in the space 
C\[-r, 0]; D(A~ 1/2 )) n C([-r, 0]; Z>(A«)), a G (|, 1). 

Proof of Theorem^. It will be shown first that (St, X) is a dissipative dynamical 
system. To that end, the below proposition is needed. 

Proposition 3. [251 Lemma 1] Let all the assumptions of Theorem [1] hold and let 
a G (§,1). Then there exists a bounded subset BV a of the space C 1 ([— r, 0]; D(A^^)) fl 
C([— r, 0]; _D(A a )), which absorbs any strong solution to the problem ([1]) and ([5]) for any 
initial function ip G C 

Second, to apply the classical theorem on the existence of a global attractor (see, for 
example (2j [3QI [6] ) , we show that (S t ,X) is asymptotically compact. Consider therefore 
any solution u(t) to the problem and (J5j) with ip G BV a as an initial function. We will 
show that for any 5 > r > and any T > 5 the set U = {u t = S t ip \ <p G BV a , t G [5, T}} 
is relatively compact in X. 

Recall that the set BV a is a ball in C\[-r, 0]; D(A~ 1/2 ))r\C([-r, 0];D(A a )) (for more 
details see [25]) and notice that, by Corollary 4 from [27], the set BV a is relatively compact 
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in C([— r, 0]; D(A 1 / 2 )) (see also [2"7l lemma 1]). It remains to show that {u(t)\ip G 
BFa, t G [5 - r, T]} is equi-continuous in C([5 -r,T); D(A~ 1/2 )). 

Proposition 4. [T9| Corollary 4.3.3 and Theorem 4.3.5]. Let A be an infinitesimal gener- 
ator of an analytic semigroup {T(t)} t >o. If / G L 1 ((0,T);F) is locally Holder continuous 
on (0,T], then for every x G Y the initial value problem 

u{t) = Au(t) + f(t),t > 0; u(0) = x 

has a unique solution u. If / G C e ([0, T]; Y"), then for every 5 > 0, Am G C e ([<5, T]; K) and 
MGC e ([5,T];F). 

Here C e ([0, T]; Y") denotes the family of all Holder continuous functions on [0, T] with the 
exponent 9 G (0, 1). In this case, Y = L 2 (fl). 

In order to apply Proposition H] to our case, we have to show that p(t) = Fi(u t ) G 
C e ([5 - r,T}; L 2 (n)) (c.f. (13^1. (1551)). Therefore, consider t G [5 - r,T] and 

| \p(t + h) - p(t) 1 1 = | |*i (u t+h ) - Fi K) 1 1 < 

< L Fl [£ BVa ] max \\A~ 1/2 {u{t + h + s) - u(t + s))\ \ < L Fl [£ BVa ] £ BVa \h\ 

s€[-r,0] 

where L^J^gyJ is the constant defined in Lemma [1] with £ BVa such that IHV'III < £ B v a 
Wip G (the existence of such £ BVa follows from Proposition [2]). Here, q = is used. 

The last inequality shows that p : [5 — r,T] — » L 2 (Q) is Lipschitz continuous, which 
is the situation to which Proposition [3] can be applied. It should also be noted that the 
family {p(t)}, for all initial ip G BV a , is uniformly Lipschitz, i.e. all the Lipschitz constants 
are lower or equal to L = Lp^BVa] • ^BV a - Then by Proposition [2J it is guaranteed (see 
the proof) that the family {u(t) \ (p G BV a , t G [5 — r, T]} is uniformly Holder continuous, 
and thus equi-continuous in C([5 — r, T]; L 2 (f2)). 

Proposition 5. [271 lemma 1] Let 5 be a Banach space. A set F of C([0, T];5) is 
relatively compact if and only if 

(i) F{t) = {fit) : / G F} is relatively compact in 5, < t < T, 

(ii) F is uniformly equicontinuous, i.e. We > 0,3r] such that \\f(t 2 ) — < e, V/ G 
F, V0 < ti < t 2 < T such that \t 2 - h\ < rj 



17 



Applying Proposition completes the proof of Theorem [2J 



□ 



As an application we can consider the diffusive Nicholson blowflies equation (see e.g. 
[29] ) with state-dependent delays, i.e. the equation (CQ) where —A is the Laplace operator 
with the Dirichlet boundary conditions, Q C R n ° is a bounded domain with a smooth 
boundary, the nonlinear (birth) function b is given by b(w) = p ■ we~ w . The function b is 
bounded, so for any delay function rj satisfying (Hl.rj), the conditions of Theorem [1] and 
Theorem |5] are satisfied. As a result, we conclude that the initial value problem ([I]) and (jSJ) 
is well-posed in X and the dynamical system (S t ,X) has a global attractor (Theorem [2]). 
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